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Abstract 



Employing profits data of Japanese companies in 2002 and 2003, we identify the non- 
Gibrat's law which holds in the middle profits region. From the law of detailed balance in all 
regions, Gibrat's law in the high region and the non-Gibrat's law in the middle region, we 
kinematically derive the profits distribution function in the high and middle range uniformly. 
The distribution function accurately fits with empirical data without any fitting parameter. 
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1 Introduction 

The distribution of wealth, income or company size is one of important issues not only in 
economics but also in econophysics [1]. In these distributions, a cumulative number N{> x) 
obeys a power-law for x which is larger than a certain threshold xq: 



This power-law and the exponent fj, are called Pareto's law and Pareto index, respectively [2] ^. 
Here x is wealth, income, profits, assets, sales, the number of employees and etc. 

The study of power-law distributions in the high range is quite significant. Because a large 
part of the total wealth, income or profits is occupied by persons or companies in the high region, 
although the number of them is a few percent. They have the possibility to influence economics. 
Power-law distributions are observed in fractal systems which have self-similarity. This means 
that there is self-similarity in economic systems. The power-law distribution in the high region 
is well investigated by using various models in econophysics. 

Recently, Fujiwara et al. [5] find that Pareto's law (and the reflection law) can be derived 
kinematically form the law of detailed balance and Gibrat's law [6] which are observed in high 
*e-mail address: ishikawa@kanazawa-gu.ac.jp 

^Recently Pareto's law is checked with high accuracy (See [3, 4] for instance). 




for X > Xq . 




region x > xq. In the proof, they assume no model and only use these two underlying laws in 
empirical data. 

The detailed balance is time-reversal symmetry: 

Pl2(xi,X2) = Pl2{x2,Xl) . (2) 

Here xi and X2 are two successive incomes, profits, assets, sales, etc, and ^12(3^1, 3^2) is a joint 
probability distribution function (pdf). Gibrat's law states that the conditional probability 
distribution of growth rate Q{R\xi) is independent of the initial value xi: 

Q{R\xi) = Q{R) . (3) 

Here growth rate R is defined as the ratio R = X2lx\ and Q{R\x\) is defined by using the pdf 
P{x\) and the joint pdf Pir{xi,R) as 

In Ref. [7], by using profits data of Japanese companies in 2002 (xi) and 2003 (2:2), it is 
confirmed that the law of detailed balance (2) holds in all regions xi > and X2 > (Fig. 1). 
We also find that Gibrat's law still holds in the extended region xi > xq and X2 > 0. Prom these 
observations, we have shown that Pareto index is also induced from the growth rate distribution 
empirically. 

On the other hand, it is also well known that the power-law is not observed below the 
threshold xq [6, 8]. Por instance, we show the profits distributions of Japanese companies in 
2002 and 2003 (Pig. 2). We find that Pareto's law holds in the high profits region and it fails 
in the middle one. The study of distributions in the middle region is as important as the study 
of power-law those. Because a large number of persons or companies is included in the middle 
region. Furthermore, it is interesting to study the breaking of fractal. 

In order to obtain the distribution in the middle region, we examine data below the threshold 
Xq. It is reported that Gibrat's law is valid only in the high region ([9] for instance). The recent 
works about the breakdown of Gibrat's law is done by Stanley's group [10]. Aoyama et al. also 
report that Gibrat's law does not hold in the middle region by using data of Japanese companies 
[11]. 

In the analysis of Gibrat's law in Ref. [7], we concentrate our attention to the region xi > xq 
and X2 > 0. In this paper, we examine the growth rate distribution in all regions xi > and 
X2 > and identify the law which is the extension of Gibrat's law ^. Wc employ profits data of 
Japanese companies in 2002 and 2003 which are available on the database "CD Eyes" published 
by Tokyo Shoko Research, Ltd. [12]. 

By using the extended Gibrat's law, we derive the distribution function in the high and 
middle profits region under the law of detailed balance. It explains empirical data with high 
accuracy. Notice that the distribution function has no fitting parameter. The parameters of the 
function are already decided in the extended Gibrat's law. 

^In this paper, we call the law in the region xi < xo and X2 > non-Gibrat's law. 
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2 Gibrat's law and non-Gibrat's law 



In this section, we reconfirm Gibrat's law in the high profits region and identify non-Gibrat's 
law in the middle one. We divide the range of xi into logarithmically equal bins as xi G 

4x [io2+0-2(n-i)^ io2+o.2n] thousand yen with n = 1, 2, • • • , 15 3. In Fig. 3, 4 and 5, the probability 
densities for r are expressed in the case of n = 1, • • • , 5, n = 6, • • • , 10 and n = 11, • • • , 15, 
respectively. The number of the companies in Fig. 3, 4 and 5 is "35, 513", "64, 205" and "25, 189", 
respectively. Here we use the log profits growth rate r = log^g R- The probability density for r 
defined by q{r\xi) is related to that for R by 

logio Q{R\xi) + r + logio(ln 10) = logio (l{r\xi) ■ (5) 

From Fig. 3, 4 and 5, we express the relation between r and g(r|a;i) as follows: 

logio = c{xi) - t+{xi) r for r>0, (6) 

logio = c{xi) +t-{xi) r for r<0. (7) 

By the use of Eq. (5) , these relations are rewritten in terms of R as 

Q{R\xi) = d{xi) for R>1 , (8) 

Q{R\xi) = d{xi) for R<1 , (9) 

with d{xi) = 10^(^1 Vln 10 . 

By applying the expressions (6) and (7) to data in Fig. 3, 4 and 5, the relation between xi 
and t±{xi) is obtained (Fig. 6). In Fig. 6, t± hardly responds to xi for n = 13,14,15. This 
means that Gibrat's law holds in the high profits region. On the other hand, t+ increases 
(decreases) linearly with log^g xi for n = 4, 5, • • • , 11. In the middle profits region, not Gibrat's 
law but the other law (non-Gibrat's law) holds as follows ^: 

t+{xi) = t+{xo) - a+ In — , (10) 

Xo 

t-{xi) = t_(xo) +a_ In— . (11) 

Xo 

In Fig. 6, a_|_ and a_ are estimated as 

q;+ ~ q;_ ~ for xi > xo , (12) 

a+ ~ a_ ~ 0.27 for Xmin < xi < Xo , (13) 

where .tq = 4 x io2+0-2(i3-i) ^ iq5 thousand yen (= 100 million yen) and a^min = 4 x 
-[^q2+o. 2(4-1) ^ 1,600 thousand yen. In this paper we call the combination of Gibrat's law 
((10), (11) and (12)) and non-Gibrat's law ((10), (11) and (13)) extended Gibrat's law. 

^In Ref. [7], we only consider the case for n = 11, 12, ■ ■ ■ , 15. 
"^In section 5, we examine another type of non-Gibrat's law. 
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3 Profits distribution in the high and middle region 

In Refs. [7], Paxeto's law (1) and the Pareto index n can be derived from the detailed balance 
(2) and Gibrat's law (3) in the high profits region. In this section, we derive profits distribution 
not only in the high profits region but also in the middle one by using the detailed balance and 

the extended Gibrat's law. 

Due to the relation of Pi2{xi, X2)dxidx2 = Pir{xi, R)dxidR under the change of variables 
from (xi,X2) to {xi,R), these two joint pdfs are related to each other, 

Pir{xi,R) =XiPi2{xi,X2). (14) 

By the use of this relation, the detailed balance (2) is rewritten in terms of Pir{xi,R) as follows: 

Pir{xi,R) = R-^Pm{x2,R-^). (15) 

Substituting the joint pdf Pir{xi, R) for the conditional probability Q{R\x\) defined in Eq. (4), 
the detailed balance is expressed as 

P{xi) ^ 1 Q{R-^\x2) , . 

P{x2) R Q{R\xi) ■ ^ ^ 

In the preceding section, the conditional probability Q{R\xi) is identified as (8) or (9). Under 
the change of variables xi X2, the conditional probability Q{R~^\x2) is expressed as 

Q{R~^\x2) = d{x2) R~^+^''^^-^ fori?<l, (17) 
Q{R^^\x2) = dix2) R^*-^''^^~^ fori?>l. (18) 

Prom the detailed balance and the extended Gibrat's law, one finds the following: 
P(xi) 



P{X2) 



= i^+*+(^l)-*-(^2) + l ^IQ-^ 

^ ^+t+(xo)-t-(a;o)+l+a+lnfi+a_ln^ ^^Q) 



for R > 1. Here we assume that the x dependence of d{x) is negligible in the derivation, the 
validity of which should be checked against the results. By expanding Eq. (20) around R = 1, 
the following differential equation is obtained 



X 

+t+{xo) — t-{xo) + 1 + (a+ + «_) In — 

•^0 



P{x) + X P'{x) = 0, (21) 



where x denotes xi. The same differential equation is obtained for < 1. The solution is given 

by 



P{x) = Cx-^'^+i) e""''' ^ (22) 

with a = («_!_ + a-) /2. Here we use the relation +i_|_(a;o) — t-{xQ) = ji which is confirmed in 
Ref. [7]. 
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4 Data fitting 



In the previous section, we derive the profits distribution function (22) in the high and 
middle profits region from the detailed balance and the extended Gibrat's law. In this section, 
we directly examine whether it fits with profits distribution data. 

In order to average the scattering of data points, we employ the cumulative number of 
companies (Fig. 2). For xi > xq, by using Eqs. (12) and (22) the cumulative distribution in the 
high profits region is expressed as 

/■oo 

Nh{> x) = N{> xo) / dt P{t) (23) 



= iV(>x„)(^-J . (24) 

On the other hand, by the use of Eqs. (13) and (22) the cumulative distribution in the middle 
profits region Xmin < xi < xq is given by 

Nm{> x) = {N{> x^in) - N{> xo)} / dt P{t) + N{> xo) (25) 

Jx 

Erf (27s) -Erf 



{N{> )-N{> xo)} 



h+2q In ^ 

fiL 

2v^ 

E^f(27s)-Erf0^^7^ 



+ N{> xo) . (26) 

Here Erf(.x) is error function defined by Erf(.T) = 2 e~^'^ d,t/ ^Jl:. 

In Fig. 7, the distribution functions (24) in the high profits region (xi > xq) and (26) in the 
middle one (xmin < a^i < a^o) explain empirical data in 2003 with high accuracy. This guarantees 
the validity of the assumption in the previous section. Notice that there is no ambiguity in 
parameter fitting, because indices ji, a and the bounds xq, Xmin is already given in the extended 
Gibrat's law ^. 



5 Another non-Gibrat's law 

In section 2, we present the non-Gibrat's law (10), (11) and (13) as a linear approximation of 
t-l-(xi) in Fig. 6. In this section, we examine another linear approximation in Fig. 8, the vertical 
axis of which is the logarithm. 

In Fig. 8, logj^Q t+ (logio increases (decreases) linearly with log^g ^1 for n = 4, 5, • • • , 11. 
This relation is expressed as 

t+(xi) = t+(xo) (27) 
t_(xi) = i_(xo) (—\ , (28) 

\XoJ 

^Pareto's law with /u = 1 is especially called Zipf's law [13]. 
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where 



t+{xo) - t^{xo) = n, (29) 
/?+ ~ 0.35 , /3_ ~ 0.36 for Xmm < xi < , (30) 
P+ ~ for xi> xo . (31) 

Here and values in section 2. From the detailed balance and this extended 

Gibrat's law, one finds 

= ^-*-(-) (^)"'-*-(-) (32) 

for R>1. 

By expanding this equation around R= 1, the following differential equation is obtained 



P(x)+a;P'(x) =0 , (33) 



where a; denotes x\. The same differential equation is obtained for i? < 1. In order to take 
— > limit into account, we rewrite this as follows: 



+M.„) ^ ( I) - i| - { (I) - lU M + 1 



P{x) + X P\x) = . (34) 



The solution is given by 

P{x) = e'^^""' ^ + 0(/?+2) + 0{pJ) (35) 

with/3 = {(3+t+{xo) + P-t-{xo)) /2. Here we use the expansion x^~-^/(a;o)^—a;~-^ = /3(lnx/a;o)/a;+ 

If we neglect 0(/3±^) terms in Eq. (35), the profits distribution function (35) can be identified 
with (22), because /3 ~ a numerically. In other words, there is no essential difference between 
two expressions (10), (11) and (27), (28) in the middle region. 



6 Conclusion 

In this paper, we have kinematically derived the profits distribution function in the high and 
middle region from the law of detailed balance and the extended Gibrat's law by employing 
profits data of Japanese companies in 2002 and 2003. 

Firstly, we have reconfirmed Gibrat's law in the high profits region {x > xq). The value of 
Xo is estimated to be about 100 million yen in Fig. 6 or 8. At the same time, we have found that 
Gibrat's law fails and another law holds in the middle profits region (xmin < x < xq). We have 
identified the non- Gibrat's Iciw cind. the value of Xmin is estimated to be about 1,600 thousand 
yen in Fig. 6 or 8. We have called the combination of Gibrat's law and the non-Gibrat's law the 
extended Gibrat's law. 
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Secondly, we have derived not only Pareto's law in the high profits region but also the 
distribution in the middle one from the detailed balance and the extended Gibrat's law. The 
derivation has been described uniformly in terms of the extended Gibrat's law. 

The profits distribution in the middle region is very similar to log-normal one: 



P{x) = Cx~^ exp 



1 1 2 ^ 



(36) 



where x is mean value and cj^ is variance. The difference between two distributions (22) and 
(36) is only the power of x. It brings a translation along the horizontal axis in the log-log plot. 
In this sense, the distribution in the middle profits region obtained in this paper is essentially 
equivalent to the log-normal one. Notice that it has no fitting parameter that the log-normal 
distribution, prepared only for data fitting, has. Indies fj,, a and the bounds Xq, Xmin is already 
given in the extended Gibrat's law. 

In the derivation of the profits distribution, we have used the detailed balance and the ex- 
tended Gibrat's law. The detailed balance is observed in a relatively stable period in economy [5]. 
The extended Gibrat's law is interpreted as follows. Companies, classified in small-scale profits 
category, have more (less) possibilities of increasing (decreasing) their profits than companies 
classified in large-scale one. In other words, it is probably difficult that companies gain large- 
scale profits in two successive years. This is the non-Gibrat's law and it leads the distribution 
in the middle profits region. 

This phenomenon is not observed above the threshold xq. For x > xq, companies, classified 
in small-scale profits category, have same possibilities of increasing (decreasing) their profits 
with companies classified in large-scale one. This is the Gibrat's law and it leads the power 
distribution in the high profits region. 

In this paper, we cannot mention the distribution in the low profits region x < .Xmin, because 
no law is observed in the region (Fig. 6 or 8). This may be caused by insufficient data in the 
low region. Lastly, we speculate possible distributions in the region. If we do not take /3± — >■ 
limit, we should solve the differential equation (33). The solution is given by 



P{x) = Cx ^ exp 



(3+ \xqJ /?_ \xoJ 



(37) 



For the case /3+ ~ 1 and (3- ~ 0, the distribution (37) takes the exponential form in Refs. [4, 14]. 
On the other hand, for the case /?+ ~ and /3_ ~ 1, the distribution (37) takes Weibull form in 
Ref. [15]. We can decide the distribution in the low profits region if suflficient data in the region 
are provided. 

We have showed that a distribution function is decided by underlying kinematics. For profits 
data we used, the distribution is power in the high region and log-normal type in the middle one. 
This docs not claim that all the distributions in the middle region arc log-normal types. For 
instance, the personal income distribution may take a different form, if the extended Gibrat's 
law changes. Even in the case, the other extended Gibrat's law will decide the other distribution 
by the use of the method in this paper. 
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Figure 1: The scatter plot of all companies, the profits of which in 2002 (xi) and 2003 {X2) 
exceeded 0, xi > and X2 > 0. The number of the companies is "132,499" . 
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Figure 2: Cumulative number distributions N(xi) and N{x2) for companies, the profits of whicfi 
in 2002 (xi) and 2003 (0:2) exceeded 0, xi > and X2 > 0. Tlie number of the companies is 
"132,499". 
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Figure 3: The probabihty density distribution q{r\xi) of the log profits growth rate r = 
logioa;2/a;i from 2002 to 2003. The data points are classified into five bins of the initial profits 
with equal mag nitude in logarithmic scale, xi G 4 x [io2+o.2(n-i) ;LQ2+o.2n] ^ 1,2, •••,5) 
thousand yen. The number of companies in this regime is "35,513" . 
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Figure 4: The probability density distribution q(r\xi) of the log profits growth rate r = 
logiQ X2/X1 from 2002 to 2003. The data points are also classified into five bins of the initial prof- 
its with equal magnitude in logarithmic scale, xi G 4 x [102+0-2(n-i) ^ ]^o^+°-2"] (n = 6, 7, • • • , 10) 
thousand yen. The number of companies in this regime is "64,205" . 
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Figure 5: The probability density distribution q{r\xi) of the log profits growth rate r = 
logiQ X2/X1 from 2002 to 2003. The data points are also classified into five bins of the initial prof- 
its with equal mag nitude in logarithmic scale, xi G 4x [lO^+^-^C"-!) , io2+0-2"] (n = 11, 12, • ■ ■ , 15) 
thousand yen. The number of companies in this regime is "25,189". 
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Figure 6: The relation between the lower bound of each bin xi and t±{xi). Prom the left, each 
data point represents n=l,2,---,15. 




Figure 7: Cumulative number distributions N{x2) for companies, the profits of which in 2002 
(xi) and 2003 (.X2) exceeded 0, xi > and X2 > 0. The cumulative distribution function derived 
from the detailed balance and the extended Gibrat's law accurately fits with the data. Indices 
/i, a and the bounds Xq, Xmin is already given in the extended Gibrat's law. 
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Figure 8: The log-log plot of the relation between the lower bound of each bin xi and t±{xi). 
Prom the left, each data point represents n = 1,2, •••,15. 
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